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Abstract 

The classical two-dimensional anisotropic triangular nearest-neighbor Ising (ATNNI) 
model is studied by the density matrix renormalization group (DMRG) technique when 
periodic boundary conditions are imposed. Applying the finite-size scaling to the DMRG 
results a commensurate-disordered (C-D) phase transition line as well as temperature and 
magnetic critical exponents are calculated. We conclude that the (C-D) phase transition 
in the ATNNI model belongs to the same universality class as the ordered-disordered 
phase transition of the Ising model. 

Analysis of semi-finite systems of small size in one or more directions has been used as 
a powerful tool in extracting of critical properties of two-dimensional classical models and 
corresponding one-dimensional quantum models. Although finite or ID systems themselves 
do not display any critical behavior, it is, however, possible to extract critical parameter values 
as well as critical exponents. Temperature, ordering magnetic field, and finite-size deviations 
from criticality are all described by the same set of the critical exponents Q. This paper 
is focused on infinite strips of finite width where the relevant numerical data are obtained 
from the transfer matrix methods, in particular, the Density Matrix Renormalization Group 
(DMRG) method. 

In 1992 the DMRG technique has been invented by S. R. White (2| in real space for one- 
dimensional (ID) quantum spin Hamiltonians. Three years later T. Nishino || applied this 
numerical technique to classical spin 2D models that is based on the renormalization group 
transformation for the transfer matrix for the open boundary conditions. DMRG treatment 
of 2D classical systems exceeds the classical Monte Carlo approach in accuracy, speed, and 
size of the systems Q| . 

Recently, we have modified the DMRG method for the 2D classical models, imposing 
periodic boundary conditions (PBC) on strip boundaries, and found a relation that helped 
to determine an optimal strip width L opt in order to obtain correct values of critical temper- 
ature and exponents H using the finite-size scaling (FSS). We have obtained results of very 
high accuracy exceeding the DMRG method with standard open boundary conditions. Our 
method does not require any extrapolation analysis of the data. 

The use of DMRG for 2D classical models may follow one of two different approaches: 

(i) DMRG method is applied to strips of finite width and from two largest transfer- 
matrix eigenvalues or the free energy estimated with high precision, the critical properties of 
the system are calculated by the FSS analysis (here, we use this approach). 



(ii) The strip width is enlarged until a steady state is reached (in the thermodynamic limit) 
when the output from the DMRG does not depend on the lattice size. Then, the DMRG yields 
properties of the 2D infinite system with spontaneously broken symmetry and mean-field-like 
behavior close to the criticality. This approach was used recently to study the high-field part 
of the ATNNI- model phase diagram ||, where approach (i) ran into convergence problems. 
We were able to show that the phase transition between the commensurate phase and the 
disordered phase proceeds via a narrow strip of an incommensurate phase. This approach 
gives also accurately the low-field part of the phase diagram, but it is not convenient for 
determination of the critical properties of the system by FSS. In distinction to the finite- 
width approach (i), it explicitly undergoes the phase transition, but its critical behavior is 
mean-field-like and the speed of calculation suffers from critical slowing-down at the phase 
transition line. Therefore, we use here approach (i) to find the low-field critical behavior of 
the ATNNI model. 

The FSS approach should give the correct critical properties of the system in the limit of 
infinite strip width. Nevertheless, it was shown in ||] that in approximate DMRG treatment 
for given size of the transfer matrix (limited by computer capacity), it is not useful to enlarge 
the strip width to too large values, because here the the DMRG results do not satisfy the 
scaling laws assumed by the FSS. Thus, an optimal width, up to which the results systemat- 
ically improve, must exist. It was also shown that the estimation of critical properties of the 
Ising and Potts models by DMRG with the periodic boundary conditions are much better 
than those with the open ones, despite the latter yields better results for the finite-width 
strips ||. 

Below the optimal strip width L opt the ratio 

R = f 2 T6(L) (1) 

is almost linear function of L while above it, it is not. {L in (1) is the width of the strip and 
Tq(L) is the critical temperature for given L.) 

The deviation of R from linearity above the optimal strip width is very fast and the ratio 
R becomes zero or infinity within enlargement of the strip by one lattice constant. Thus, if 
R = or R — > oo (i.e. if the numerator or the denominator tends to zero or changes its sign), 
we accept that L as the strip width for further calculations and call it the optimal width L opt 
of the strip. The critical temperature for the optimal width Tfj(L opt ) is taken as the best 
approximation of the critical temperature of the 2D system studied, and at this temperature 
the critical exponents of the system are calculated. The values of the critical exponents are 
sensitive to Tq and must be determined with the due care. 

In the FSS approach, the critical exponents are derived from the scaling behavior of the 
correlation length and free energy at critical point, where they depend on strip width L in 
the following way [jij: 

K\~ L < ] Kl-Vti' CL ~L<'- d (2) 

where and K\ are the derivatives of inverse correlation length K with respect to tem- 
perature T and second derivative with respect to ordering (magnetic) field h, respectively, 
and cl is the specific heat, i.e. the second derivative of the free energy with respect to tem- 
perature. The two temperature exponents y^ and should be equal to each other. The 
exponents and yh determine the critical behavior of all statistical quantities characterizing 
the system. The critical exponents of specific heat, magnetization and correlation length can 
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be calculated from t/t and yt as follows: a = 2 — (3 = 2 y ^ h , v = y^ 1 ■ Other critical expo- 
nents can be obtained from the scaling equations y& = f3 + 7, 7 = (3{5 — 1), and rj = 2 — 72/T 

Further, we demonstrate the capabilities of our approach to find the critical properties of 
2D spin lattice model on Ising model with different symmetries of the lattice, where critical 
temperatures and critical indices are known from exact solutions, and ATNNI model where 
the phase diagram is generally unknown and the critical indices are predicted from symmetry 
considerations. 

The 2D classical anisotropic triangular nearest-neighbor Ising (ATNNI) model is given by 
the Hamiltonian 

n = - J I Yl a i a i+s + aa i a i+3 J - H J2 °i ( 3 ) 

1 \&=l,2 J l 

with the antiferromagnetic coupling J < and spins U{ = ±1. The numbers 1, 2, and 3 
are lattice directions in the ATNNI model. The coupling J is multiplied by the parameter a 
(0 < a < 1) along the direction 3 as depicted in Figure [l] (a). 

ATNNI model on the triangular lattice 
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Figure 1: (a) The triangular lattice of the ATNNI model, (b) The phase diagram of the 
ATNNI model for a=04 obtained by DMRG g]. 

This model was studied by Domany and Schaub [H and in ||, and it was shown that 
its phase diagram, as a plot of temperature T and external magnetic field H (for a = 0.4), 
exhibits four different phases: two commensurate phases (I) and {IT}, a disordered phase, 
and an incommensurate phase, see Figure [l] (b). Commensurate phase (J) occurs at mag- 
netic field H < 2.4. This structure satisfies the Lifshitz condition, and it is characterized by 
a one-dimensional representation of the lattice symmetry group, i.e. its phase transition is 
predicted to belong to the Ising universality class §. Domany and Schaub tried to confirm 
this prediction by numerical calculation of the exponent yx> but due to the low-order approx- 
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imation it differed from the expected value by more than 10% and the magnetic exponent 
was not calculated at all. 

We have calculated critical properties of the ATNNI model at the phase transition line be- 
tween the commensurate (I) and disordered phase. To illustrate the accuracy of the method, 
we have calculated critical properties of the exactly solvable models: ferromagnetic and an- 
tiferromagnetic Ising model on square and triangular lattices at zero magnetic field as well 
as the zero-magnetic-field ATNNI model which critical temperature is given by the equation 

We have used the FSS analysis of DMRG results with superblock consisting of 8 Ising spins 
and 4 multi-spin variable, each acquiring 85 values (m = 85). The computational effort at 
this approximation is less than for the classical transfer matrix method of strip width equal 
to 17 lattice constants. However, the DMRG enables to treat wider strip (of tens of lattice 
constants) up to the optimal width further improving the values of the critical parameters. 

The first, important step of the calculations is determination of the critical temperature 
Tq, see Table EL of which the best estimation for given m is T<S(L opt ) calculated from FSS 
approach ||. At this temperature the values of the critical exponents are derived from the 
scaling laws (2). 



Table 1: Critical temperatures T<S obtained from (1) with DMRG compared with the exact 



ones T, 



(exact) 



C 



The symbols □ and A describe square and triangular lattices, respectively. 



model 


H 


J c 


rp (exact) 


□ Ising 


0.0 


2.2691851 


2.2691853 


□ AF Ising 


0.0 


2.2691848 


2.2691853 


A Ising 


0.0 


3.640955 


3.640957 


A ATNNI 


0.0 


1.55352 


1.55362 


A ATNNI 


0.5 


1.52867 


unknown 


A ATNNI 


1.0 


1.45135 


unknown 


A ATNNI 


1.5 


1.31105 


unknown 


A ATNNI 


2.0 


1.07009 


unknown 



As the quantities appearing in (2) are first and second derivatives of the free energy and 
correlation length, the effect of approximation starts to manifest at lower strip width than 
L opt . The criterion determining strip width at which the value of the critical exponent may 
be still acceptable, was taken completely analogous to that for critical temperature, Eq. (1). 
The accepted values of the critical exponents are denoted by filled symbols in Figs || (a) and 
(b). 

The critical exponent yx is determined more precisely from the free energy y^f than from 

the correlation length y™ \ as for the evaluation of the former only the largest eigenvalue of 
the superblock matrix is needed in distinction to the correlation length, to whose calculation 
the ratio of the largest and the second largest eigenvalue is necessary. This point is irrelevant 
for the models with a symmetric transfer matrix (Ising models in Table [l]), but significant for 



the ATNNI model with a non-symmetric transfer matrix [11]. The plot of thermal critical 
exponent vs. strip width is shown in Figure § (a). For increasing lattice size they both 
tend to the Ising value 1. The convergence also depends on the magnetic field. It gets worse 
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(a) (b) 

Figure 2: (a) The plot of thermal critical exponents y^ for different fields in the ATNNI 

model, (b ) The plot of magnetic critical exponents . The filled symbols denote the accepted 
critical exponents satisfying Eq. (1) 

for magnetic field close to the multi-critical point H = 2.4. Here the reliability of the DMRG 
breaks down at rather small strip width, as well. The accepted values depicted by black 
symbols are listed in Table |2[ 



Table 2: Critical exponents of various 2D spin models calculated by the DMRG method 
with PBC and FSS analysis. The exact critical exponents of the Ising models are as follows: 
2/t = 1 and y^ = 1.875. 



model 


H 


y T 






a 




V 


□ Ising 


0.0 


1.0000009 


0.99999994 


1.875002 


0.0000017 


i 

8.00012 


1.00000006 


□ AF Ising 


0.0 


1.0000009 


0.99999994 


1.875126 


0.0000017 


8.00804 


1.00000006 


A Ising 


0.0 


1.0000014 


0.99999943 


1.875030 


0.0000027 


8.00192 


1.00000057 


A ATNNI 


0.0 


1.0000022 


0.9947 


1.87005 


0.000004 


7^0 


1.00527 


A ATNNI 


0.5 


1.0000280 


0.9902 


1.87098 


0.000056 


TIE 


1.00993 


A ATNNI 


1.0 


1.0000580 


0.9902 


1.87062 


0.000116 


7~73 


1.00993 


A ATNNI 


1.5 


1.0000767 


0.9911 


1.86939 


0.000153 


THE 


1.00893 


A ATNNI 


2.0 


0.9998366 


1.0122 


1.86902 


0.000327 


7153 


0.98795 



The critical exponent y^ describes the decay of the order parameter at the phase tran- 
sition line from the commensurate phase (I) to the disordered phase. The structure (/} 
consists of two ferromagnetically ordered sublattices each with different magnetization. As 
the external magnetic field H is generally non-zero in ATNNI model, the total magnetization 
(sum of both sublattice magnetizations) is non-zero, as well. The difference between the two 
magnetization is taken as the order parameter in this case. The small ordering field h used 
for calculation of the derivative acquires opposite sign at each of the two sublattices. The 
accuracy of the calculations of the magnetic exponent is smaller than that of the thermal 
exponent in case of exactly solvable models listed in Table ^. Thus, we can expect a lower 
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accuracy also for ATNNI model. All the exponents depicted in Fig. |2| (b) are below 1.871. 
Extrapolations to L — > oo for H = 0.5-1.5 give values of about 1.872, i.e. (5 = y^j, 
which still differs from the Ising value y^ = 1.875 and corresponding j3 = g. Note that the 
value of y^ is extremely sensitive to the correct determination of the critical temperature. 
A very small decrease of its value would shift y^ to the expected Ising value. At modest 
magnetic field, where our calculations are assumed to be more accurate, the plots of y^ for 
different magnetic field lie on the same curve what suggests that not only i s a universal 
quantity independent of H, but the corrections to it for finite L are universal, as well. 

In conclusion, it can be stated that the DMRG method with periodic boundary conditions 
reproduces with a high accuracy the critical properties of exactly solvable models and confirms 
the prediction that the C-D phase transition for magnetic fields H = 0-2.4 belongs to the 
universality class of the Ising model. 

A. G. thanks to T. Nishino for the fruitful correspondence and help. Most of this work 
has been done by Compaq Fortran on the HPC Alpha 21264 in Kobe University in Japan. 
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